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Abstract 

In multielectron bubbles, the electrons form an effectively two-dimensional layer at the inner 
surface of the bubble in helium. The modes of oscillation of the bubble surface (the ripplons) are 
influenced by the charge redistribution of the electrons along the surface. The dispersion relation 
for these charge redistribution modes ('longitudinal plasmons') is derived and the coupling of these 
modes to the ripplons is analysed. We find that the ripplon-plasmon coupling in a multielectron 
bubble differs markedly from that of electrons a flat helium surface. An equation is presented 
relating the spherical harmonic components of the charge redistribution to those of the shape 
deformation of the bubble. 

PACS numbers: 73.20.-r, 64.70.Dv, 68.35. Ja, 47.55.Dz 



1 



I. INTRODUCTION 



Multielectron bubbles are 0.1 /iin - 100 /xm sized cavities inside liquid helium, that contain 
helium vapor at vapor pressure and a nanometer-thick electron layer anchored to the surface 
of the bubble jH- The dependence of the radius of the bubble on the pressure applied to 
the bubble allows to vary the surface density of the spherical two-dimensional electron 
system in the bubble over nearly four orders of magnitude. As such, these objects offer 
the prospect to study the two-dimensional electron system on helium in regimes hitherto 
inaccessible. In particular, this would allow to study the two-dimensional electron Wigner 
crystal at densities not currently achievable for electrons on a flat helium surface. 

Early experimental evidence for 2D Wigner crystallization of electrons on a liquid-He 
surface Q relied on the detection of the coupled plasmon-ripplon modes 0|. InRef. 0,the 
ripplon modes characteristic for a bubble inside liquid helium were calculated and investi- 
gated as a function of the pressure applied on the bubble. In this communication we extend 
this investigation and take into account the redistribution of charge along the bubble surface 
when the bubble deforms. This allows us to derive the dispersion relation for the spherical 
'plasmon' modes and the plasmon-ripplon coupling peculiar for the spherical surface. 



II. SMALL AMPLITUDE DEFORMATIONS AND CHARGE REDISTRIBUTION 

The surface of a bubble in liquid helium can be described by a function R{Q) which gives 
the distance of the helium surface from the center of the bubble, in a direction given by the 
spherical angle Q = {d,(f)}. In equilibrium (at zero pressure) the multielectron bubble has 
a spherical surface with radius i?b- We expand the deformation away from the spherical 
surface in spherical harmonics: 

oo e 

R{Q) = Rh + J2J2 QimYe„^{n). (1) 

i=l m=-l 

where i?b is the angle-averaged radius of the bubble, Qnm is the amplitude of the spherical 
harmonic deformation mode indexed by {£, m}, and Y^m,(^) is the corresponding spherical 
harmonic. These modes of deformation are referred to as ripplons in analogy with the 
surface modes of a fiat surface of liquid helium. In Ref. Q the dispersion relation for the 
ripplons in a multielectron bubble has been studied as a function of pressure and related to 
the bubble stability. In this communication,we extend the results of in order to take the 
coupling and mixing of the ripplons and the electronic modes into account. 

In a multielectron bubble, the electrons are confined to a thin (1 nm) layer at the inner 
surface of the bubble. This layer is anchored to the surface of the bubble, so that when the 
surface deforms, the layer conforms to the new bubble surface. The electrons can redistribute 
themselves inside the spherical layer, so that the surface density of electrons is no longer 
uniform. We describe the surface density with the function ns{^l) giving the number of 
electrons in a solid angle dfl around Q : 

N °° ^ 

""^^^^ " 4^ + 5Z nimYlm{^). (2) 

b i=i m=-l 
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This function is normalized such that N = J ns{fl)RldQ is the total number of electrons in 
the bubble or on the droplet. The n^m's represent the amplitudes of charge redistributions 
corresponding to spherical harmonics Y^m- 

The total potential energy of the multielectron bubble or droplet can be separated into 
several contributions: (i) a term from the surface tension energy Us = crS, with a ~ 
3.6 X 10~^ J/m^ and S the surface of the deformed bubble or droplet; (ii) a pressure related 
term Up = pV, with p the difference in pressure outside and inside the bubble and V the 
volume of the bubble; (iii) a term representing the electrostatic energy Uc of the electron 
layer. In this last term, the quantum corrections (such as the exchange contribution) to the 
electrostatic energy can be neglected 0. The first two terms, Us and Up, have been derived 
ini: 

oo e 

Us = AnaRl + ^ ^ ^ (^^ + £ + 2) \Q,^f , (3) 
£=1 m=-e 

A OD £ 

Up = -^pRl+pRbY. Yl l^^-l'- (4) 

i=l m=-i 

Whereas Us is an exact expression. Up is an expansion up to second order in Qim and 
the third order term has been neglected. In what follows, we assume small amplitude 
deformations and small amplitude charge redistributions such that 



VW+^\Qem\ « Rb, (5) 



VWTT)M « (6) 



The electrostatic potential V^(r) of the deformed MEB with a non-uniform surface electron 
density is calculated straightforwardly within the framework of the Maxwell equations by 
expanding the potential inside and outside the deformed bubble in their respective spherical 
decompositions and imposing the electrostatic boundary conditions at the surface. The 
potential energy associated with this electrostatic potential is given by 

Uc = -{1/2) j ns{n)v{R,n)dn. (7) 

If we keep only terms up to second order in both n^m and Qem, we find 



i=l m=-£ 

^ \^ \ ^ —62(1+ 1) I „ |2 



=1 m=-£ 

00 



=1 m=-£ 



where e is the electron charge, ei is the dielectic constant of the medium inside the surface, 
and £2 is the dielectric constant of the medium outside the surface. For a multielectron 
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bubble, 62 = e with e = 1.0572 the dielectric constant of liquid helium, and ei = 1 (if the 
vapour pressure of the helium is low enough, the bubble has vacuum inside). 

The ripplon contribution to the kinetic energy of the MEB is associated with the motion 
of the liquid helium surface. Following a derivation of Lord Rayleigh for oscillating liquid 
droplets p, we find 



2 ^ ^ 

1=1 m=-i 



Q 



(9) 



where p ~ 145 kg/m^ is the density of liquid helium. Note that for a bubble, i + 1 appears 
in the denominator instead of i for the droplet in Lord Rayleigh's treatment. 
The 'classical' kinetic energy of the electrons is given by 

N 

T, = !^5^w^(r,) (10) 

i=i 

with nie the electron mass, and w(rj) represents the displacement of electron j out of its 
equilibrium position r^. Since we assumed we can use the formula 



T 



ill 



I w\r)dh (11) 

J surface 



and express the field of displacements w(r) as a sum of a longitudinal field WL(r) and 
transverse field WT(r). We investigate the effect of the longitudinal field, which can be 
written as a gradient of a scalar potential. The (divergence-free) transverse field is not 
considered here since it does not couple to the ripplons. Using 

V-WL(r) = l-^n<,(^,0) (12) 
we can express the kinetic energy of the electrons as 

1 47rmei?g , . ,2 

^e = 2LL]w(£TT)'"^'"' ^''^ 

e=i m=-l ^ ' 

We have checked that our approach yielded, for the 3D and the 2D electron gas, the known 
expressions for the plasma frequencies. 



III. RESULTS AND DISCUSSION 

The full Lagrangian of the bubble, including the ripplon modes and charge redistribution 
modes, is given by substituting expressions ©jllini),®,®,® in £ = Tj. + Te — Us~Up — Uc- 
The result can be brought in the following form 



+ 5Z (^^™ ~ ^li^Hm) + ^^^ernQim- (14) 
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This is the central result of this communication. We will now proceed to discuss the result 
term by term. The first term in expression (fT^ contains the Lagrangian describing the 
radial motion (the anharmonic breathing mode): 

Cr = f^Rl - -pRt - AnaRl - (15) 

In what follows we assume that Rb = and the bubble radius is given by its equilibrium 
value. 



A. Ripplon modes 

The Lagrangian C also contains a part representing the harmonic oscillation of the ripplon 
modes, with oscillator mass term 

(16) 



and bare ripplon frequency: 



+ 1- 




a{& + i+l)+pRb- 



2„2 



eii + 1] 



AneRl 



(17) 



The bare ripplon frequencies and their dependence on the pressure were the subject of jlj. 
Note that, in the limit of very large bubbles, R^ ^ oo, the dispersion relation for ripplons 
on a fiat helium surface |7] is recovered: if the momentum k is identified with i/Rb, then 
the typical k^^'^ is retrieved. Note furthermore that, if Rb is set equal to the equilibrium 
coulomb radius, the modes £=1,2 obtain a zero frequency in agreement with the result of 
Salomaa and Williams For typical multielectron bubbles with radii in the range of 1-100 
/xm (and numbers of electrons of the order 10^ — 10^) the ripplon frequencies for i < 1000 
lie typically in the MHz-GHz range. 



B. Longitudinal plasmon modes 

Next, the Lagrangean ()14j) contains a part representing the (small amplitude) harmonic 
oscillation of the classical charge redistribution around its equilibrium. The oscillator mass 
term for these oscillations is 

and the corresponding frequencies are given by 



In the limit of large bubbles, this frequency corresponds to that of the longitudinal plasmon 
frequency of a 2D electron system on the helium surface If the electrons form a Wigner 
lattice, this is the longitudinal phonon frequency. If the momentum k is identified with i/Rb, 
then the typical k^^'^ disperson for this longitudinal plasmon is retrieved. Based on this 
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correspondence, we will refer to these modes in a multielectron bubble as the '(longitudinal) 
plasmon modes of the MEB' (or 'longitudinal phonon modes of the MEB' if the electrons 
form a Wigner lattice). 

For typical multielectron bubbles with radii in the range of 1-100 /im (and numbers of 
electrons of the order 10'^ — 10^) the longitudinal plasmon frequencies for £ < 1000 lie typically 
in the GHz-THz range. For bubbles with N > 10^, the longitudinal plasmon frequencies are 
larger than the ripplon frequencies, u!p{£) ^ u!r{£)- 

The equations of motion for the electron charge redistribution on a deformed bubble are 
given by: 

d dC _ dC 

dt dhim driim 

4^ mghim = -meUp{i)nem + liQim- (20) 

When ijjpijl) ^ the electrons can redistribute much faster on the surface of the bubble 
than the bubble surface can deform. Thus, we can make an adiabatic approximation to 
find the charge redistribution of the electrons on a bubble with a given deformation. We 
find from ()20|) that if the bubble surface deformation is described by a given set of spherical 
deformation amplitudes Qim-, the equilibrium charge distribution on the deformed bubble 
must satisfy the relation 

nim (^+l)Q^m 



with rio = A^/(47ri?^). 



C. Ripplon - longitudinal plasmon coupling 

The last term in the Lagrangean (|Tl|l represents the coupling between the longitudinal 
plasmons and the ripplons, with coupling strength given by 

For electrons on a flat helium surface, such a coiipling between ripplons and longitudinal 
plasmons/phonons was derived by Fisher et al. j4| and detected experimentally ]|3|. For 
electrons on the inner surface of a deformed bubble, we find a similar coupling, but only 
ripplon and longitudinal plasmon modes which have the same angular momentum couple to 
each other. After the diagonalisation of the ripplon-plasmon part of the Lagrangian ([njl . 
we arrive at the eigenfrequencies, 

^^1,2 (^) = Hi^) + ± \\/Hif^) - ^m' + 47|. (23) 

In Fig. 1, we show the eigenfrequencies fli {£) and Q2 {£) for N = 10^ as a function of £ 
and p. The frequency Q2 {£) is close to the ripplon frequency derived for an MEB within the 
approximation of a conducting surface Consequently, this frequency can be attributed 
to the ripplon modes, renormalized due to the ripplon-plasmon mixing. The admixture with 
the longitudinal plasmon mode is weak. The other branch of oscillations with the frequencies 
fl2 (0 t)e related to the longitudinal plasmon mode admixed with a small component 
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of ripplon nature. In typical multielectron bubbles, the difference in frequency between 
the ripplon and the plasmon modes weakens the coupling between these modes. Unlike for 
electrons on a fiat surface Q the mixing is weak. Finally, note that the multielectron bubble 
will be stable when both flf and fl"^ positive. This condition is equivalent to 

u;l{i)Lu',{i) > 7£ (24) 

By substituting the results (fT^ . (fTTj) and in the inequality it can easily be seen 
that when the radius of the bubble is equal to the equilibrium radius, the ripplon-plasmon 
mixing does not change the criterion of stability for a MEB formulated in Ref. 



IV. CONCLUSIONS 

In this communication we studied the longitudinal plasmon-ripplon modes in a multi- 
electron bubble. The central result is the Lagrangean (fT^ with the ripplon and plasmon 
frequencies fll7|) . (fT^ and the coupling strength (j^^ . Unlike for electrons on a fiat surface, we 
find that the mixing of the modes is weak for typical multielectron bubbles, because the bare 
longitudinal plasmon and ripplon frequencies are different {ujp{i) 3> ujr{i))- The conditions 
for bubble stability discussed in Ref. Q remain unaltered by the ripplon-plasmon mixing. 
The analytic expressions for the longitudinal plasmon frequencies and the coupling strength 
as a function of the angular momentum allow us to derive a compact formula describing 




Fig- 1 



FIG. 1: Vibrational eigenfrequencies of an MEB in liquid helium, fii (l) (phonon modes) and 
^^2 (0 (renormalized ripplon modes) given by Eq. ((2S1) for N = 10^ as a function of / and of p. 
Inset: a schematic picture of the directions of motion for phonons and ripplons in an MEB. The 
ripplons are excitations of the helium surface (with typical frequencies in the MHz-GHz range), 
while the phonons are related to the motion of electrons (with typical frequencies in the THz 
regime) tangential to the bubble surface. Note the different scales for Qi and 1^2. 
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the charge redistribution taking place on a deformed bubble. For a bubble deformation in 
a particular spherical harmonic mode {£,m}, the charge will redistribute itself according to 
the same spherical harmonic mode, and the amplitudes of deformation will be proportional 
to each other. 
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